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ABSTRACT 

Feedback from the central black hole in active galactic nuclei (AGN) may be respon- 
sible for establishing the observed Mbh — cr relation and limiting the bulge stellar 
mass of the host galaxy. Here we explore the possibility of AGN feedback triggering 
star formation in the host galaxy. We consider a shell of dusty gas, driven outwards 
by radiation pressure, and analyse its escape/trapping condition in the galactic halo 
for different underlying dark matter potentials. In the isothermal potential, we obtain 
that the standard condition setting the observed Mbh — o relation is not sufficient to 
clear gas out of the entire galaxy; whereas the same condition is formally sufficient in 
the case of the Hcrnquist and Navarro-Frcnk- White profiles. The squeezing and com- 
pression of the inhomogeneous interstellar medium during the ejection process can 
trigger star formation within the feedback-driven shell. We estimate the resulting star 
formation rate and total additional stellar mass. In this picture, new stars are formed 
at increasingly larger radii and successively populate the outer regions of the host 
galaxy. This characteristic pattern may be compared with the observed 'inside-out' 
growth of massive galaxies. 

Key words: black hole physics - galaxies: active - galaxies: evolution - stars: forma- 
tion 



1 INTRODUCTION 

Black holes in active galactic nuclei (AGN) respond to the 
accretion process by feeding back energy to the surroundings 
in both radiative and kinetic forms. The interaction between 
AGN feedback and the ambient material can have a great 
impact on scales of the host galaxy and beyond. The energy 
released by accretion onto the central black hole clearly ex- 
ceeds the binding energy of the host galaxy bulge, and it 
is now widely agreed that feedback from the central AGN 
plays an important role in the formation and evolution of 
the host galaxy. However, the question is how this energy 
and momentum resulting from the accretion process is cou- 
pled and interacts with the surrounding medium. This can 
in principle be achieved in a variety of ways, involving both 
mechanical and radiative processes, such as jets, winds, and 
radiation pressure. 

Observations indicate a number of direct relationships 
between the mass of the central black hole and the char- 
acteristic properties of the host galaxy. One of the tight- 
est relations is given by the empirical correlation observed 
between the black hole mass and the bulge stellar velocity 
dispersion, of the form Mbh oc a 4 , known as the Mbh — o~ 
relation (Gebhardt et al. 2000; Ferrarese & Merritt 2000). A 
similar relation holds between the black hole mass and the 
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bulge stellar mass, leading to Mbh ~ 10 3 M b ui g e (Marconi 
& Hunt 2003; Haring & Rix 2004). The existence of such 
correlations suggests a close coupling between the evolution 
of the central black hole and that of its host galaxy. 

The observational correlations can be interpreted in 
terms of AGN feedback, and different feedback models have 
been proposed and discussed in the literature (Silk & Rees 
1998; Fabian 1999; Fabian et al. 2002; King 2003, 2005; 
Murray et al. 2005; Fabian et al. 2006, 2008; King 2010; 
Silk & Nusser 2010; King et al. 2011). A simple deriva- 
tion of the Mbh — o relation is obtained by considering 
momentum balance, which leads to a scaling of the form 
M B h = Jq2^ o" 4 (= M CT ), where / g is the gas mass fraction 
and ctt is the Thomson cross section (e.g. Fabian 2012, and 
references therein) . This is in remarkable agreement with the 
observational relation. It is generally assumed that once the 
central black hole reaches the critical mass, M a , the result- 
ing feedback blows away the surrounding material from the 
galaxy bulge, eventually clearing gas out of the host galaxy. 
The removal of the ambient material, which forms the po- 
tential fuel for the feeding of the central black hole and the 
gas reservoir for star formation, prevents further accretion 
and formation of stars. This leads to a switch-off of the black 
hole growth and a quenching of star formation in the host 
galaxy, as seen in some numerical simulations of galaxy evo- 
lution (e.g. Springel et al. 2005; Di Matteo et al. 2005). Thus 
the black hole seems to self-regulate its own growth, and at 
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the same time limit the stellar mass content of the host by 
expelling surrounding material. According to this standard 
view, the connection between the central black hole and its 
host galaxy is essentially via a negative feedback mechanism. 

However, feedback from the central black hole might 
play additional roles in the evolution of its host galaxy. Here 
we explore the possibility of AGN feedback triggering star 
formation in the host galaxy (Fabian 2012). The interstellar 
medium of a galaxy is composed by cold and dense molecular 
gas, which form potential sites for star formation. Radiation 
from the central black hole is expected to be absorbed by 
dust embedded in the gas. We consider a shell of dusty gas 
being swept outwards by radiation pressure, and follow its 
evolution in the galactic potential. We assume that all of 
the radiation momentum output is available for driving the 
shell, and we ignore the details of radiative transfer. The 
squeezing and compression of the inhomogeneous interstel- 
lar medium may trigger star formation within the outflow- 
ing shell. Stars that are formed will eventually be dropped 
out of the expanding shell, and start to follow their own or- 
bits in the gravitational potential of the host galaxy. Several 
such episodes of star formation, due to AGN feedback ac- 
tivity, may be envisaged during the lifetime of a galaxy. At 
each episode, mass is deposited at outer radii with a fraction 
of gas being converted into stars. This process adds stellar 
mass in the outer regions of the galaxy, thus contributing 
to the growth of the host. In order to test the plausibility 
of such a mechanism, we start by examining more in detail 
the behaviour of the feedback-driven shell in different grav- 
itational potentials due to the contributions of dark matter 
and old stars. 

The paper is organized as follows. We first discuss dif- 
ferent physical mechanisms possibly driving AGN feedback 
(Sect. 2). We next examine the escape/trapping conditions 
of the outflowing shell in different underlying dark matter 
potentials (Sect. 3). In Sect. 4 we explicitly introduce star 
formation in the expanding shell, and estimate the star for- 
mation rate along with the corresponding additional stellar 
mass. The resulting implications are discussed in Sect. 5. 



2 WINDS AND RADIATION PRESSURE 

As mentioned in the Introduction, AGN feedback can be 
driven by different physical mechanisms. One way of driv- 
ing large-scale outflows from the central nucleus is via high- 
velocity winds (King 2003, 2005). Indeed there is observa- 
tional evidence of winds and outflows detected in a number 
of objects (Pounds et al. 2003; Sturm et al. 2011; Tombesi 
et al. 2012). The wind is assumed to be launched from the 
immediate vicinity of the central black hole, with an outflow 
momentum of the same order of the Eddington limited pho- 
ton momentum: M w v ~ Le/c. As the wind is generated in 
the small central region, the mass flow rate is limited and the 
outflow velocity must be very high. Based on the assumption 
that the wind mass flow rate is equal to the Eddington mass 
accretion rate (M w ~ Me), one obtains that v/c ~ 77, where 
77 ~ 0.1 is the standard radiative efficiency (King 2010; 
Zubovas & King 2012). The derived high velocity (v ~ 0.1c) 
corresponds to the escape speed at a launch radius of a few 
hundred gravitational radii. When such high-velocity winds 
encounter the cold interstellar medium of the host galaxy, 



strong shocks will inevitably result. The hot shocked gas is 
cooled via inverse Compton emission; but the cooling pro- 
cess is only efficient in the inner regions, where the cool- 
ing time is shorter than the flow time. Beyond a critical 
radius (of the order of ~ 1 kiloparsec) there is no way to ef- 
ficiently cool the shocked gas, and the initially momentum- 
driven flow, considered further by McQuillin & McLaughlin 
(2012), becomes energy-driven. The subsequent evolution of 
the energy-driven outflow is extensively analysed in King 
(2010) and King et al. (2011). 

We do not further pursue this issue here, and next con- 
sider the direct effects of radiation pressure on dusty gas, 
as first introduced for AGN feedback by Fabian (1999) and 
discussed in Fabian et al. (2002). In this case, the main in- 
teraction is attributed to absorption of the central radiation 
by dust grains embedded in the gas. The radiative momen- 
tum is transferred to the gas through the coupling between 
dust and gas, which is further enhanced in the central en- 
vironment close to the ionising source where both dust and 
gas are partially charged (Murray et al. 2005; Fabian et al. 
2008). This strong coupling ensures that gas can be effi- 
ciently dragged out with the dust. We note that sublimation 
of dust is relevant only within the inner ~parsec region (on 
scales smaller than the broad line region of the central AGN, 
cf Murray et al. 2005) and should not be an issue on the 
scales discussed here. We assume that the dusty gas which is 
being driven by radiation pressure is strongly radiating and 
remains cold (T <C 10 4 K), while outer gas swept up by the 
expanding shell is shocked to high temperatures (~ 10 6 — 10 7 
K) and subsequently cools down. Provided that a significant 
amount of the dust survives sputtering and other destruc- 
tive processes in the shocks (for a discussion on dust grain 
survival, see Jones & Nuth 2011), the mass of swept up ma- 
terial in the shell steadily increases. A critical radius can 
be defined, beyond which the outflowing shell becomes opti- 
cally thin to the central radiation. The exact location of the 
critical radius depends on the coupling process, i.e. electron 
scattering or dust absorption. It is well known that the dust 
cross section (ad) is much larger than the Thomson cross 
section (ar), and the ratio is found to be of the order of 
(Jd/cfT ~ 1000 (Fabian 2012, and references therein). Dust 
opacities are typically of several hundred cm 2 /g (Murray 
et al. 2005), and the corresponding critical radius is of the 
order of a few tens of kiloparsecs. Thus radiation pressure 
on dust grains is another viable process of driving gas out 
to large radii, and in the following we consider the dynamics 
of flows driven by radiation pressure on dusty gas. 



3 ESCAPE/TRAPPING CONDITIONS 

We assume that radiation pressure on dusty gas sweeps up 
the surrounding material into a shell, and follow its tem- 
poral evolution in the gravitational potential of the host 
galaxy. We start by examining the conditions for the escape 
(or trapping) of the outflowing shell from the galactic halo, 
which is a prerequisite for the further study of star forma- 
tion triggered within the feedback-driven shell. The general 
form of the equation of motion of the shell is given by: 

at c H 
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where L is the luminosity of the central source, M s (r) the 
enclosed gas mass, and Mdm(V) the dark matter mass. We 
assume that at sufficiently large radii (r 3> GA/bh/o" 2 ), the 
black hole contribution can be neglected such that the po- 
tential is dominated by the dark matter component. In the 
following, we consider the evolution of the expanding shell in 
three types of underlying dark matter potentials, and derive 
the shell escape/trapping conditions in a direct and simple 
scheme. A more complex and detailed analysis of the dynam- 
ics of momentum-driven shells can be found in McQuillin & 
McLaughlin (2012). 



3.1 Isothermal potential 

Modelling the galaxy as an isothermal sphere, the density 
profile of dark matter is given by 



p(r) 



(2) 



2nGr 2 

where a is the velocity dispersion. The enclosed mass inside 
a sphere of radius r is 

(3) 



/2a 2 
p(r)r 2 dr = -^-r , 



The gas mass is assumed to be a fraction f s of the dark 
matter mass 



Af g (r) = / g M DM (r) 



2f e a 2 
G 



(4) 



In the isothermal case, an analytical solution to the equation 
of motion can be derived. Integrating twice we obtain: 



r (t) = \ r 2 + 2r v t + 



GL 
2/g«r 2 



- 2a 2 t 2 



(5) 



where ro and Vq are the initial position and velocity, respec- 
tively. We note that this solution is analogous to the result 
previously obtained by King (2010) for the shock pattern 
far from the central black hole. The velocity is given by 
v(t) — ^y. The velocity of the shell can also be expressed as 
a function of radius as: 



v(r) 



GL 



2a 2 



1 r o 



(6) 



The escape or trapping condition of the shell in the 
galactic potential is mainly determined by a critical lumi- 
nosity. The critical luminosity is defined by considering the 
equilibrium between outward radiation pressure and inward 
gravitational forces, which leads to 



4/ g ca 4 



G 



(7) 



We see that in the isothermal case the critical luminosity is 
constant and independent of radius. The numerical value is 
of the order of L c = 4.6 x 10 46 / g ,o.i60" 4 oo er g s_1 f° r typical 
values of the velocity dispersion, a — 0200 x 200 km s , and 
gas mass fraction, / g = / g ,o.i6 x 0.16. For luminosities lower 
than the critical value (L < L c ), the outward force due to 
radiation pressure is not sufficient to compensate the inward 
gravitational attraction: the shell falls back, independently 
of the initial conditions. But once the central source exceeds 
the critical luminosity (L ^ L c ) the shell is able to expand 
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Figure 1. Evolution of the shell in the isothermal potential for 
different luminosities, close to the critical value: L ~ 0.9 L c (black 
dotted line), L ~ 0.95 L c (blue dashed line), L ~ 1.01 L c (green 
solid line). The initial conditions are the same for the three curves: 
ro = 100 pc, vo = 200 km s — 1 . The shell can expand outwards 
only for L ^ L c . 



outwards indefinitely. The evolution of the shell for differ- 
ent central luminosities is shown in Fig. 1. In Fig. 2 we plot 
the radial velocity profile of the expanding shell for different 
initial conditions. We see that at large radii (r ^> ro) the 
expanding shell always reaches a constant velocity, indepen- 
dently of the initial conditions. From Eq. (6) we can derive 
the value of the asymptotic velocity, Vao = v(r 2> ro): 



GL 
V^o~ 2 



2a 2 



(8) 



We observe that the asymptotic velocity has no dependence 
on the initial conditions, and is primarily determined by 
the central luminosity. The shell velocity is larger for higher 
luminosities, scaling as i>oo oc yL. The velocity of the out- 
flowing shell can be compared with the local escape velocity. 
For an isothermal sphere, the escape velocity is constant and 
independent of radius: v csc — 2a. The formal condition for 
the escape of the shell is given by 



GL 
Zhca 2 



-2a 2 \ >2a. 



(9) 



The escape condition (Eq. 9) then implies that the central 
luminosity should exceed a critical value: 



L>L' 



12/ g «T 4 
G 



3L C 



(10) 



This defines another critical luminosity, which is equal to 
three times the critical luminosity for the shell expansion 
derived in Eq. (7). It implies that the luminosity of the cen- 
tral source should at least reach L' c — 3L C for the shell to 
be able to escape the galactic halo. 



3.1.1 Implications for the Mbh 



relation 



The luminosity of the central source is related to the mass 
of the central black hole. If one assumes the standard 
Eddington limit for accretion onto the black hole 
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Figure 2. Radial velocity profile of the expanding shell in the 
isothermal potential. The central luminosity is L ~ 9.4 X 10 46 erg 
s _1 (~ 2 L c ); the three curves indicate different initial conditions: 
ro = 100 pc, vo = 500 km s — 1 (black dashed line), rrj = 10 pc, 
vg = 200 km s -1 (green dash-dot line), ro = 50 pc, vq = km 
s" 1 (blue dotted line). At large radii the shell reaches a constant 
asymptotic velocity, independent of the initial conditions. 



Figure 3. Radial velocity profile of the expanding shell (black 
solid line) and local escape speed (cyan dashed line) in a Hernquist 
potential. The escape velocity is defined as d csc = ^/2|$(r)|, 
where the gravitational potential is obtained by integrating the 
Poisson equation V 2 4> = AnGp. The central luminosity is L c ~ 
4.6 X 10 46 erg s _1 . The shell velocity exceeds the local escape 
velocity at very large radii. 
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where n — ox/m-p is the electron scattering opacity, the 
escape condition (Eq. 9) implies: 



M BH > = 3M CT 



(12) 



Note that k in the above expression refers only to the inter- 
action of radiation and matter limiting accretion onto the 
black hole, and not to the interaction with the interstellar 
medium of the host galaxy. We thus recover the standard 
Mbh — o relation, but with an additional factor of three 
compared to the usual value of M a ■ This suggests that the 
condition Mbh = M CT only allows a shell to expand towards 
larger radii, but is not sufficient for the shell to be expelled 
from the host galaxy. In order for the shell to completely es- 
cape the galaxy halo, the black hole mass must exceed 3M CT 
(see also McQuillin & McLaughlin (2012)). 

3.2 Hernquist profile 

Galaxies are embedded in dark matter haloes, which dom- 
inate the outer gravitational potential. Dark matter haloes 
are not ideal isothermal spheres and are better approxi- 
mated by profiles with varying slopes. We first consider the 
Hernquist density profile (Hernquist 1990) given by 

Mdm r a 



Pdm(»~) 



(13) 



2-7T r(r + r a ) 3 ' 

where Mdm is the total mass, and r a is a characteristic scale 
length where the density profile changes shape. Compared 
to the isothermal case, the Hernquist density profile is shal- 
lower at small radii (p oc 1/r) and steeper in the outer re- 
gions (p oc 1/r 4 ). The corresponding mass distribution is 
given by 



M DM (r) = Mdm 



(r + r a ) 2 



(14) 



The gas mass is taken to be a fraction / g of the dark matter 
mass. By analogy with the isothermal case, a critical lumi- 
nosity can be defined by equating the radiation pressure and 
the gravitational force terms: 



L c ,h = £ c (r a ) = c/ { 



CM DM 
16r| 



(15) 



As in the previously discussed isothermal case, once the cen- 
tral luminosity exceeds the critical value, shells can expand 
outwards to arbitrarily large radii. In a Hernquist potential, 
the shell velocity reaches a minimum around the scale ra- 
dius and increases at larger radii, eventually exceeding the 
local escape velocity (Fig. 3). The shape of the radial veloc- 
ity profile implies that the shell is able to escape the galactic 
halo, once the luminosity reaches the critical value given in 
Eq. (15). Thus in a Hernquist potential, the critical luminos- 
ity required for the shell escape is equal to that defined for 
the shell expansion. However, in order to escape the galactic 
halo, the shell needs to reach very large radii where its veloc- 
ity starts to increase and eventually exceed the local escape 
velocity. This implies that AGN feedback should be power- 
ful enough to bring the shell at the required large distances. 

In the above discussion we have implicitly assumed that 
the cold gas directly traces the dark matter distribution. But 
this assumption may not be adequate, in particular in the 
inner regions. The stellar density distribution generally fol- 
lows p oc 1/r 2 , and cold gas from which stars ultimately 
form should also follow a similar distribution. Therefore one 
may separate the two contributions in the equation of mo- 
tion, and assume an isothermal distribution for the gas along 
with a Hernquist profile for the dark matter component. The 
resulting velocity profile is shown in Fig. 4. We observe that 
the shell does not accelerate to very high velocities at small 
radii, and is similar to the isothermal case in the inner re- 
gions; while at outer radii the velocity profile follows the 



AGN feedback and triggering of star formation 5 




in 

radius[kpc] 



E 

4- 750 



in 

radius[kpc] 



Figure 4. Radial velocity profile of the expanding shell (black 
solid line) and local escape speed (cyan dashed line) in a poten- 
tial composed by isothermal gas and Hernquist dark matter. The 
central luminosity is L c ~ 4.6 X 10 46 erg s — 1 . The shell velocity 
exceeds the local escape velocity at very large radii. 



Figure 5. Radial velocity profile of the expanding shell (black 
solid line) and local escape speed (cyan dashed line) in a NFW 
potential. The central luminosity is L c ~ 1.5 X 10 46 erg s — 1 . The 
shell velocity exceeds the local escape velocity at very large radii. 



characteristic shape of the Hernquist potential, which po- 
tentially allows the shell escape. 



3.3 NFW profile 

We next consider the so-called NFW profile (Navarro, Frenk 
& White 1996, 1997), assumed to be a universal profile 
describing dark matter haloes. The density distribution is 
given by 

pcSc 



p(r) = 



(r/r s )(l + r/r s ) 2 



(16) 



where p c = is the critical density, S c is a characteristic 
density, and r s is a scale radius. At large radii (r ;§> r s ), 
the NFW profile is shallower (p oc f/r 3 ) than the Hernquist 
profile; while both profiles tend to p oc l/r in the central 
regions. The associated mass profile is given by 



M(r) = 4iv5 c p c ri 



In 1 + 



r 

r + r s 



(17) 



We repeat the same analysis as for the Hernquist poten- 
tial case, and first determine the critical luminosity defined 
as the maximum of the L c (r) function. As in the previously 
discussed cases, expanding shells can reach large radii only 
for luminosities exceeding the critical value. The velocity 
profile of the NFW potential is qualitatively similar to the 
Hernquist case: the shell velocity reaches a minimum beyond 
which it starts to increase, eventually exceeding the local 
escape velocity at very large radii (Fig. 5). Considering an 
isothermal gas distribution and a NFW dark matter gives a 
qualitatively similar result to the corresponding case in the 
Hernquist potential. 



4 TRIGGERED STAR FORMATION 

In the previous sections we have analysed the evolution of 
the radiation pressure driven shell in different gravitational 



potentials. The surrounding medium is swept up by the ex- 
panding shell; the squeezing and compression of the cold gas 
induced by the passage of the shell may cause local density 
enhancements. This could, in turn, lead to a triggering of 
star formation within the outflowing shell. We adopt a sim- 
ple prescription for the triggered star formation in the shell 
and estimate the star formation rate as: 



Mi, 



M s (r) 
tflow(r) 



(18) 



where is the star formation efficiency, and ta OVJ (r) = 
is the local flow time. The observed star formation efficiency 
is typically of the order of a few percent, and somewhat 
higher in starburst systems. In the isothermal case, the gas 
mass scales with radius and Eq. (18) leads to 



Mi, 



- J %-«(r) ■ 



(19) 



We observe that the star formation rate directly scales with 
the velocity of the shell, Mi, oc v(r). We have seen that in the 
isothermal case the shell velocity tends towards a constant 
asymptotic value at large radii. This implies that the star 
formation rate is also constant at large radii: 



M* i00 = M*(r -> oo) ~ e, 



2/gg 2 
G 



(20) 



As «oo oc yL oc %/M, more luminous or equivalently more 
massive objects tend to have higher star formation rates. For 
typical galaxy parameters and shell velocities, we obtain star 
formation rates of the order of M* )O0 ~ 10 — 100 Mq yr , 
for ~ 0.01 — 0.1, respectively (Fig. 6). Once formed, stars 
decouple from the expanding shell; the initial velocity of the 
newly formed stars corresponds to the shell velocity at the 
formation radius. As long as the initial speed is lower than 
the local escape velocity, stars that are formed will remain 
bound to the galaxy. Thus star formation in the outflow- 
ing shell contributes to the increase of the host stellar mass 
in the outer regions. Integrating over the total AGN feed- 
back lifetime, one may obtain a rough order-of-magnitude 
estimate for the total additional stellar mass formed in this 
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Figure 6. Star formation rate as a function of radius in the 
isothermal potential. The central luminosity is L ~ 9.2 X 10 46 
erg s — 1 (~ 2 L c ); the two curves indicate different star forma- 
tion efficiencies: = 0.01 (black solid line), e* = 0.1 (red dashed 
line). The star formation rate is constant at large radii. 
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Figure 8. Star formation rate as a function of radius for an 
isothermal gas distribution and Hernquist dark matter profile. 
The line symbols are the same as in Fig. 6. 



isothermal gas and NFW dark matter is again similar to the 
corresponding case for the Hernquist potential. 
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Figure 7. Star formation rate as a function of radius in the 
Hernquist potential. The central luminosity is L ~ 9.2 X 10 46 
erg s (~ 2 L c ), and the line symbols are the same as in Fig. 6. 



way. In the isothermal case, the total stellar mass can be of 
the order of several 10 9 Mq for one AGN episode lasting a 
Salpeter time (~ 5 x 10 7 yr). 

The same computations can be performed in the case of 
the Hernquist potential. The resulting star formation rate is 
shown in Fig. (7); the associated stellar mass is comparable 
to the isothermal case. For the combined case of isothermal 
gas and Hernquist dark matter, the radial profile of the star 
formation rate in the inner regions is very similar to the 
isothermal case (Fig. 8), and the total additional stellar mass 
is also of the same order. The star formation rate in the 
case of a NFW potential seems to be somewhat lower, with 
correspondingly lower stellar masses. The combined case of 



5 DISCUSSION 

It is usually assumed that when the central black hole 
reaches the critical mass, M a , the resulting feedback blows 
away the surrounding material, shutting off further black 
hole growth and suppressing star formation. However, we 
have seen that this condition is not sufficient to completely 
expel gas from the host galaxy in the case of an isothermal 
potential. In fact, for Mbh = M CT , the outflowing shell ve- 
locity tends to a constant asymptotic value at large radii, 
which is lower than the local escape velocity. The shell thus 
remains trapped in the galactic halo, and once the AGN 
activity turns off, it may eventually fall back (providing po- 
tential fuel for further accretion and star formation). The 
actual condition for the shell escape, i.e. Voo ^ 2a, requires 
the black hole mass to exceed Mbh = This is three 

times the critical value of the standard Mbh — a relation. 
It has been argued that black holes lying on the observa- 
tional A/bh — & relation cannot expel material from the 
galactic halo by radiation pressure alone, and that an ad- 
ditional momentum boost is required (Silk & Nusser 2010). 
The shortfall seems to be a factor of a few, which is consis- 
tent with our result. Recent numerical simulations of radia- 
tion pressure-driven feedback also show that for L ~ L c , gas 
is only evacuated from the central nuclear region, while the 
critical luminosity needs to be exceeded by a factor of sev- 
eral in order to clear gas out of the entire galaxy (Debuhr 
et al. 2010, 2011). However, if the galaxy is residing in a 
group or cluster environment, gas in the outer regions could 
be removed by ram pressure or tidal stripping, and thus en- 
vironmental effects also play an important role. On the other 
hand, if an initially momentum-driven wind switches to an 
energy-driven form, then a black hole mass exceeding the 
Mbh — o value may not be required (King et al. 2011). 

Previous studies of AGN feedback generally assumed 
an isothermal profile for both the gas and the dark mat- 
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ter distributions (Silk & Rees 1998; Fabian 1999; Fabian 
ct al. 2002; King 2003, 2005; Silk & Nusser 2010). Recently, 
momentum-driven feedback in non-isothermal dark matter 
haloes has also been analysed (McQuillin & McLaughlin 
2012). Here we follow the evolution of the feedback-driven 
shell and determine its escape/trapping conditions in differ- 
ent underlying dark matter potentials. At large radii, the 
Hernquist and NFW density profiles are steeper than the 
isothermal one, and the resulting radial velocity profiles are 
also qualitatively different: the shell velocity increases and 
eventually exceeds the local escape velocity at very large 
radii. Therefore escape is formally possible for M = M CT 
in a Hernquist or NFW potential, provided that the shell 
is able to reach such great distances. But this in turn re- 
quires extremely long timescales for the feedback process, 
and one should take into account the AGN lifetime. For in- 
stance in the case of a Hernquist potential, for M = M a 
the shell would take ~ 10 9 yr to reach the critical radius 
(where its velocity becomes larger than the local escape ve- 
locity) , which is very long compared to typical AGN activity 
time scales. For higher central luminosities, the critical ra- 
dius may be reached on more plausible time scales of the 
order of a few Salpeter times (~ 10 s yr) . 

Feedback from the central black hole is usually invoked 
in galaxy formation and evolution models to suppress star 
formation. However AGN feedback may also operate in the 
opposite direction and even trigger star formation in the 
host galaxy. In fact, triggering of star formation associated 
with radio jet activity, has been proposed in the past to 
explain the observed alignment of radio and optical struc- 
tures in high-redshift radio galaxies (Rees 1989; Begelman & 
Cioffi 1989). A positive feedback mechanism, induced by the 
AGN jet and leading to enhanced star formation, has also 
been considered as a source powering luminous starbursts 
(Silk 2005). Numerical simulations of radio jet-induced star 
formation have also been recently performed (Gaibler et al. 
2012). 

Here we explicitly consider triggered star formation in 
the radiation pressure-driven shell, where the characteris- 
tic time scale is given by the local flow time fe ow (r). We 
recall that we are considering cold, dusty, and primarily 
molecular gas from which stars are expected to form. As 
the shell moves on, stars are gradually dropped out of the 
expanding shell, and subsequently follow their own orbits de- 
termined by the gravitational potential of the galaxy. The 
newly formed stars are expected to follow nearly radial or- 
bits, although the exact stellar orbit depends on the ini- 
tial conditions, i.e. formation radius and local shell veloc- 
ity. Observations of bright elliptical galaxies show that the 
velocity dispersion profile steeply declines towards the ex- 
terior (Romanowsky et al. 2003; Gerhard 2010, and refer- 
ences therein) . This indicates that the orbits of the stars in 
the outer halos of bright ellipticals are dominated by radial 
motion, which may be consistent with our interpretation. 
Overall, star formation in the outflowing shell suggests that 
new stars are formed at increasingly larger radii, successively 
populating the outer regions of the host galaxy. 

This characteristic pattern may be compared with ob- 
servational studies of the outer growth of galaxies. Detailed 
observations of the growth of massive galaxies from redshift 
z ~ 2 to the present have shown that the total stellar mass 
of these galaxies has increased by a factor of ~2 over this 



time range (van Dokkum et al. 2010). But this doubling does 
not occur in a homogeneous way throughout the galaxy: the 
stellar mass in the central regions remains roughly constant, 
while the increase in mass mainly occurs at outer radii. This 
suggests that the mass growth takes place in the outer re- 
gions of the galaxy with a gradual build-up of the stellar 
mass at larger radii, implying a sort of 'inside-out' growth. 
Comparison of a sample of compact, quiescent galaxies at 
z ~ 2.3 with nearby massive elliptical galaxies suggests that 
the compact high-redshift objects may form the central cores 
of today's ellipticals (Bezanson et al. 2009). Recent mea- 
surements also confirm the extremely small sizes of massive 
quiescent galaxies at high redshift (Szomoru et al. 2012). 
All these observations seem to support the inside-out growth 
scenario, and point toward an evolutionary scenario whereby 
the galaxy growth mainly takes place in the outer regions, 
possibly by accretion of stars through merger processes. 

Alternatively, the observed outer growth of massive 
galaxies may be interpreted in terms of AGN feedback, and 
the above sketched picture of star formation in the feedback- 
driven shell seems at least qualitatively consistent with ob- 
servations. However, we should mention that observational 
studies usually analyse the quiescent phases of galaxies in se- 
lected samples at given redshifts; whereas we discuss active 
phases when the central AGN is switched on and is trig- 
gering star formation. Observationally such activity phases 
may be better classified as quasars or ultra-luminous in- 
frared galaxies. One can consider several episodes of AGN 
activity and associated feedback. Accretion of material, via 
some form of cold flows, should replenish the gas content of 
the host galaxy. Gas added to the galaxy will then fuel both 
the central black hole and feedback-driven star formation; 
the black hole mass increases and so does the stellar veloc- 
ity dispersion. As a result of the repetition of many such 
episodes, galaxies will stay on the Mbh — a relation. 

It is also interesting to note that recent observations 
indicate the existence of massive galaxies at z ~ 2 which 
are much more compact compared to galaxies of compara- 
ble mass at z ~ 0. These findings imply a significant evo- 
lution in size, with an increase of a factor of ~ 5 in effec- 
tive radius, but without a correspondingly strong increase in 
mass. Different physical interpretations of the size evolution 
of galaxies have been considered, including major mergers, 
minor mergers, and adiabatic expansion (Naab et al. 2009; 
Bezanson et al. 2009; Fan et al. 2008). The observed size 
evolution may also be interpreted in the framework of the 
AGN feedback-driven star formation. The next step will be 
to analyse more in detail the fate of the newly formed stars 
by computing their subsequent orbits in the gravitational 
potential of the host galaxy. This should allow us to make 
more quantitative comparisons with the observed 'inside- 
out' growth of massive galaxies. A succession of AGN feed- 
back episodes should also be considered and its potential 
role in the build-up and evolution of the host galaxy. 
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